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These days, as high energy particle colliders become unavailable for testing speculative theoret¬ 
ical ideas, physicists are looking to other environments that may provide extreme conditions where 
theory confronts physical reality. One such circumstance may arise at high temperature T, which 
perhaps can be attained in heavy ion collisions or in astrophysical settings. It is natural therefore 
to examine the high-temperature behavior of the standard model, and here I shall report on recent 
progress in constructing the high-T limit of QCD. 

In studying a held theory at hnite temperature, the simplest approach is the so-called imaginary- 
time formalism. We continue time to the imaginary interval [0, l/zT] and consider bosonic 
(fermionic) helds to be periodic (anti-periodic) on that interval. Perturbative calculations are 
performed by the usual Feynman rules as at zero temperature, except that in the conjugate energy- 
momentum, Fourier-transformed space, the energy variable (conjugate to the periodic time 
variable) becomes discrete — it is 2'nnT, (n integer) for bosons. From this one immediately sees 
that at high temperature — in the limiting case, at inhnite temperature — the time direction 
disappears, because the temporal interval shrinks to zero. Only zero-energy processes survive, 
since “non-vanishing energy” necessarily means high energy owing to the discreteness of the en¬ 
ergy variable p^ ~ 27rnT, and therefore all modes with n # 0 decouple at large T. In this way 
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a Euclidean three-dimensional field theory becomes effective at high temperatures and describes 
essentially static processes. [|| 

Let me repeat in greater detail. Finite-T, imaginary-time perturbation theory makes use of 
conventional diagrammatic analysis in “momentum” space, with modified “energy” variables, as 
indicated above. Specifically a spinless boson propagator is 



D{p) - pQ- in{2n)T 

Pq — p — 

(1) 

while 

a spin-i fermion propagator reads 



S{p) - po - tTT{2n + l)T 

— '~f • p — m 

(2) 

The 

zero-temperature integration measure / becomes 

replaced by 


( 2 n)^ • Thus it is seen that Bose exchange between two 0(g) vertices contributes 

'En=-oo f 1^9 S where g is the coupling strength. In the large T limit, all n / 0 

terms (formally) vanish and only the n = 0 term survives. One is left with f 9 ^^qVT. 

This is a Bose exchange graph in a Euclidian 3-dimensional theory, with effective coupling g^/T. 
Similar reasoning leads to the conclusion that fermions decouple at large T. 

While all this is quick and simple, it may be physically inadequate. First of all, frequently one is 
interested in non-static processes in real time, so complicated analytic continuation from imaginary 
time needs to be made before passing to the high-T limit, which in imaginary time describes only 
static processes. Also one may wish to study amplitudes where the real external energy is neither 
large nor zero, even though virtual internal energies are high. 

Another reason that the above may be inadequate emerges when we consider massless fields 
(such as those that occur in QCD). We have seen that the n = 0 mode leaves a propagator that 
behaves as ^ when mass vanishes, and a phase space of d?p. It is well known that this kind of 

kinematics at low momenta leads to infrared divergences in perturbation theory even for off-mass- 
shell amplitudes — Green’s functions in massless Bosonic field theories possess infrared divergences 
in naive perturbation theory. Since physical QCD does not suffer from off-mass-shell infrared 
divergences, perturbation theory must be resummed. 

Thus the formal arguments for the emergence of a 3-dimensional theory at high-T need be 
re-examined for QCD. Nevertheless, even if unreliable, the arguments alert us to the possibility 
that 3-dimensional field theoretic structures may emerge in the high-T regime. Indeed this occurs, 
although not in a direct, straightforward fashion; this will be demonstrated presently. 

Here is a graphical argument to the same end discussed above: viz. the need to resum pertur¬ 
bation theory. Consider a one-loop amplitude Hi {p ), 

ni(p) = j dk h{p,k) , 
given by the graph in the figure. 
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Compare this to a two-loop amplitude 112 (p), 

^2{p) = J dkl2{p,k) , 

in which Hi is an insertion, as in the figure below. 


n 2 (p) 



Following PisarskijQ we estimate the relative importance of 112 to IIi by the ratio of their inte¬ 
grands, 


n2 h 
Ih^h 


^ A:2 


Here g is the coupling constant, and the k"^ in the denominator reflects the fact that we are 
considering a massless particle, as in QCD. Clearly the —> 0 limit is relevant to the question 
whether the higher order graph can be neglected relative to the lower order one. Because one finds 
that for small k and large T, ni(/i;) behaves as T^, the ratio n 2 /ni is j I?. As a result when k 
is 0{gT) or smaller the two-loop amplitude is not negligible compared to the one-loop amplitude. 
Thus graphs with “soft” external momenta [0{gT) or smaller] have to be included as insertions in 
higher order calculations. 

A terminology has arisen: graphs with generic/soft external moment [0{gT) where g is small 
and T is large] and large internal momenta [the internal momenta are integration variables in an 
amplitude; when T is large they are 0(T), hence also large] are called “hard thermal loops.” P, Q] 
Much study has been expended on them and finally a general picture has emerged. Before presenting 
general results, let us look at a specific example — a 2-point Green’s function. 

It needs to be appreciated that in the imaginary-time formalism the correlation functions are 
unique and definite. But passage to real time, requires continuing from the integer-valued “energy” 
to a continuous variable, and this cannot be performed uniquely. This reflects the fact that in real 
time there exists a variety of correlation functions: time ordered products, retarded commutators, 
advanced commutators, etc. Essentially one is seeing the consequence of the fact that a Euclidean 
Laplacian possesses a unique inverse, whereas giving an inverse for the Minkowskian d’Alembertian 
requires specifying temporal boundary conditions, and a variety of answers can be gotten with a 
variety of boundary conditions. 

Thus, when presenting results one needs to specify precisely what one is computing. 

We shall consider a correlation function for two fermionic currents, in the 1-loop approximation. 




-i{f{x)f{y)) 

r d^k 


(2vr) 




{k) 


The QCD result differs from the QED result by a group theoretical multiplicative factor, so we 
present high-T results only for the latter, in real-time, and consider the time-ordained product 
as well as the retarded commutator . 
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possesses a real and an imaginary part. It is found that at large T, the real parts of 
and coincide. 


rp2 T'^h'^ 

-Ren'^'^(A:) = + 


6 *2 

where the projection operators are 


\k^ 


1 H—ttt 111 

2\k\ 


- \k\ 


fco + lfcl 


1 

3 ^ 2 2 


= g^^^-k^k^jk^ 
0 


PH'^ = 


if /r or z/ = 0 
(5*1 — PW l\k^\ otherwise 


For the imaginary part, which is present only for space-like arguments, different expressions 
found 


are 


—Imn^'^(A:) = —'Kff’^{k) 
k^ 




.2^ ^ _|_ 


-lmPi>^^{k) = 


irk"^ 


r-0 

T 


3 ^ 2 

e{-k'^ 


-pr + -pr 

3 ^ 2 ^ 


A unified presentation of these formulas is achieved in a dispersive representation. For the 
retarded function this reads 

while the time-ordered expression is 

U'^(k) = + J dk'o + gfco/T _ fc) 

The dispersive expressions may also be used to give the imaginary-time formula. 


Y[i^^ (]A — iji*** ik?\ -I- f dk' 

11 imaginary I ^ J y _ 2T:inT 

In all the above formulas, is a real subtraction term. 

Note that a universal statement about high-T behavior can be made only for the absorptive part 
: it is 0{T‘^). This also characterizes but possesses an additional 0{T^) imaginary 
part, which is seen to arise from the additional term in 11^^^ involving the bosonic distribution 
function ^kQ/T_i ■ Finally, the Il^maginary amplitude has a temperature behavior determined by its 

time 

external “energy” = 2'KinT. If this is replaced by a fixed ko (T-independent) or if only the n = 0 
mode is considered, then one may assign an 0(T'^) behavior to this quantity as well. 

In conclusion, we assert that the 2-point correlation function behaves as O(T^), where it is 
understood that this statement is to be applied to the retarded amplitude, or to the imaginary 
time amplitude with its “energy” argument continued away from 27rmr.|||] 

Similar analysis has been performed on the higher-point functions and this work has culminated 
with the discovery (Braaten, Pisarski, Frenkel, Taylor) |^, of a remarkable simplicity in their 

structure. To describe this simplicity, we do not discuss the individual n-point functions, but 
rather their sum multiplied by powers of the vector potential, viz. we consider the generating 
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functional for single-particle irreducible Green’s functions with gange field external lines in the 
hard thermal limit. (Effectively, we are dealing with continued imaginary-time amplitndes.) We 
call this qnantity EhtlC^) and it is compnted in an SU{N) gauge theory containing Np fermion 
species of the fundamental representation. Ehtl is found (i) to be proportional to (A^-|- ^Np), (ii) 
to behave as at high temperatnre, and (iii) to be gauge invariant. 


A WvQvW A -|- 




A A 



+ ' 


= iN + ^Np)^TnTLiA) 


TnTLiU-^AU + U-^dU) =rHTL(^) 


r-l 


(Henceforth g, the coupling constant, is scaled to unity.) A fnrther kinematical simplification in 
Ehtl has also been established. To explain this we define two light-like four-vectors depending 
on a unit three-vector q, pointing in an arbitrary direction. 


q-q = l , Q±Q±i, = 0 , = 1 

Coordinates and potentials are projected onto Q^. 

= Xf.Q'^ , d± = Q±-^ > A± = Af.Q'^ 

The additional fact that is now known is that (iv) after separating an ultralocal contribution from 
ThtL; the remainder may be written as an average over the angles of g of a functional W that 
depends only on A_|_; also this functional is non-local only on the two-dimensional x^ plane, and 
is ultralocal in the remaining directions, perpendicnlar to the x^ plane. [“Ultralocal” means that 
any potentially non-local kernel k{x,y) is in fact a J-fnnction of the difference k{x,y) oc 5{x — y).] 

rHTL(A) = 2pJ d^x A^{x)A^{x) + J dQgW{A+) 


These results are established in perturbation theory, and a pertnrbative expansion of VU(A_|_), i.e. a 
power series in A_|_, exhibits the above mentioned properties. A natnral question is whether one 
can sum the series to obtain an expression for W{Ap). 

Important progress on this problem was made when it was observed (Taylor, Wong)|^ that 
the gauge-invariance condition can be imposed infinitesimally, wherenpon it leads to a functional 
differential equation for IU(A+), which is best presented as 


d 5 
dx+ 


W{A+) + ljd^x A^+(x)A^(x) 


d 

—dx~ 


[A“ ] + 



W(A+) + 1 


d^x A^(x)A^(x) 


= 0 
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In other words we seek a quantity, call it 


S{A+) = i I d^xAl{x)Al{x) , 

which is a functional on a two-dimensional manifold {x"*",®”}, depends on a single functional 
variable and satisfies 

“1” = x+ , “2” = -X- , A1 = Al 

Another suggestive version of the above is gotten by defining A 2 = j^- 

diA^ - d2A1 + f^^^A\Al = 0 

To solve the functional equation and produce an expression for W (A+), we now turn to a completely 
different corner of physics, and that is Chern-Simons theory at zero temperature. 

The Chern-Simons term is a peculiar gauge theoretic topological structure that can be con¬ 
structed in odd dimensions, and here we consider it in 3-dimensional space-time. 

/cs oc J (fx Tr {d^ApA^ -I- ^A^AfsA^) 

This object was introduced into physics over a decade ago, and since that time it has been put 
to various physical and mathematical uses. Indeed one of our originally stated motivations for 
studying the Chern-Simons term was its possible relevance to high-temperature gauge theory. Q 
Here following Efraty and Nair,[^ we shall employ the Chern-Simons term for a determination of 
the hard thermal loop generating functional, Thtl- 

Since it is the space-time integral of a density, /cs may be viewed as the action for a quantum 
field theory in (2-|-l)-dimensional space-time, and the corresponding Lagrangian would then be 
given by a two-dimensional, spatial integral of a Lagrange density. 

/cs oc J dt Lcs 

Lqs oc j d?x (^A^Ai AqFi 2 ^ 

I have separated the temporal index (0) from the two spatial ones (1,2) and have indicated time 
differentiation by an over dot. F 12 is the non-Abelian field strength, defined on a two-dimensional 
plane. 

F^2 = - asA? + r^^A{A^2 

Examining the Lagrangian, we see that it has the form 

L pq- XH{p,q) 

where Af plays the role of p, Af that of q, F 12 is like a Hamiltonian and Aq acts like the Lagrange 
multiplier A, which forces the Hamiltonian to vanish; here Aq enforces the vanishing of Fi 2 - 

Fi2 = 0 

The analogy instructs us how the Chern-Simons theory should be quantized. 
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We postulate equal-time commutation relations, like those between p and q. 

A?(r), 4(r')] =iS^^6{r-r') 

In order to satisfy the condition enforced by the Lagrange multiplier, we demand that F 12 , operating 
on “allowed” states, annihilate them. 

F,^2\ ) = 0 

This equation can be explicitly presented in a Schrodinger-like representation for the Chern- 
Simons quantum field theory, where the state is a functional of Af. The action of the operators Af 
and A 2 is by multiplication and functional differentiation, respectively. 

I ) ~ 

TJI ) ~ 

' i5A<l ^ 

This, of course, is just the field theoretic analog of the quantum mechanical situation where states 
are functions of q, the q operator acts by multiplication, and the p operator by differentiation. In 
the Schrodinger representation, the condition that states be annihilated by F^2 

{diA<i-d2A\ + UcA\A'2) I ) = 0 
leads to a functional differential equation. 

If we dehne S' by T = we get equivalently 

This equation comprises the entire content of Chern-Simons quantum field theory. S is the Chern- 
Simons eikonal, which gives the exact wave functional owing to the simple dynamics of the theory. 
Also the above eikonal equation is recognized to be precisely the equation for the hard thermal 
loop generating functional, given above. 

Let me elaborate on the connection with eikonal-WKB ideas. Let us recall that in particle 
quantum mechanics, when the wave function ^{q) is written in eikonal form 

^{q) = F^ig) 

then the WKB approximation to S{q) is given by the integral of the canonical 1-form pdq 

S{q) = J p{q)dq 

where p{q), the momentum, is taken to be function of the coordinate q, by virtue of satisfying the 
equation of motion. 

^ + = E 

p{q) = J 2 E- 2V{q) 
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Analogously, in the present field theory application, the eikonal 5'(Ai) may be written as 

5(Ai)= A^^iA[)VAf 

where A^iAi) is functional of Ai determined by the equation of motion 

diA^ - d2Af + r'’^A\Al = 0 

Since, by construction = A 2 , it is clear that as a consequence S satisfies the required equation. 

However, we reiterate that here there is no WKB approximation: everything is exact. 

The gained advantage for thermal physics is that “acceptable” Chern-Simons states, i.e. solu¬ 
tions to the above functional equations, were constructed long ago,|9| and one can now take over 
those results to the hard thermal loop problem. One knows from the Chern-Simons work that T 
and S are given by a 2-dimensional fermionic determinant, i.e. by the Polyakov-Wiegman expres¬ 
sion. While these are not described by very explicit formulas, many properties are understood, and 
the hope is that one can use these properties to obtain further information about high-temperature 
QCD processes. 

For example one can give a very explicit series expansion for Fhtl in terms of powers of A 

Thtl = i / rW,,.4.4 + i I T>Sl^AAA + ■■■ 

(i) 

where the non-local kernels F^^p^ nre known explicitly. This power series may be used to systematize 
the resummation procedure for pertubative theory. Here is what one does: perturbation theory 
for Green’s functions may be organized with the help of a functional integral, where the integrand 
contains (among other factors) where /qcd is the QCD action. We now rewrite that as 

^*|^QCD(^)+^rHTL(^)-^rHTL(^)| 


where m = T Obviously nothing has changed, because we have merely added and sub¬ 
tracted the hard-thermal-loop generating functional. Next we introduce a loop counting parameter 
1: in an /-expansion, different powers of I correspond to different numbers of loops, but at the end 
I is set to unity. The resummed action is then taken to be 


g^Tesummed —— g f 






One readily verifies that an expansion in powers of / describes the resummed perturbation theory, 
free of infrared divergences. 

Even though the closed form for Fhtl is not very explicit, a much more explicit formula can 
be gotten for its functional derivative ^^5^- This may be identified with an induced current, 

which is then used as a source in the Yang-Mills equation. Thereby one obtains a non-Abelian 
generalization of the Kubo equation, which governs the response of the hot quark gluon plasma to 
external disturbances. ^ 

2 

TYl 

r) 

2 -/induced 

From the known properties of the fermionic determinant — hard thermal loop generating functional 
— one can show that jj^duced given by 

induced = / ^ {Q+ («-(^) - + Q- | 








where a± are solutions to the equations 

5+a_ — A+ + [A+, a_] = 0 

d+A- — d-a+ + [a+, A^] = 0 

Evidently Jinduced! determined by the above equations, is a non-local and non-linear functional 
of the vector potential A^. 

There now have appeared several alternative derivations of the Kubo equation. Blaizot and 
Iancu|T^ have analyzed the Schwinger-Dyson equations in the hard thermal regime; they truncated 
them at the 1-loop level, made further kinematical approximations that are justified in the hard 
thermal limit, and they too arrived at the Kubo equation. Equivalently the argument may be 
presented succinctly in the language of the composite effective action, which is truncated at 
the 1-loop (semi-classical) level — two-particle irreducible graphs are omitted. The stationarity 
condition on the 1-loop action is the gauge invariance constraint on Thtl- Einally, there is one 
more, entirely different derivation — which perhaps is the most interesting because it relies on 
classical physics. |]l2| We shall give the argument presently, but first we discuss solutions for the 
Kubo equation. 

To solve the Kubo equation, one must determine a± for arbitrary A±, thereby obtaining an 
expression for the induced current, as a functional of A±. Since the functional is non-local and 
non-linear, it does not appear possible to construct it explicitly in all generality. However, special 
cases can be readily handled. 

In the Abelian case, everything commutes and linearizes. One can determine a± in terms of 

A±. 

a± = 

Incidentally, this formula exemplifies the kinematical simplicity, mentioned above, of hard thermal 
loops: the nonlocality of l/d± lies entirely in the plane. With the above form for a± 

inserted into the Kubo equation, the solution can be constructed explicitly. It coincides with the 
results obtained by Silin long ago, on the basis of the Boltzmann-Vlasov equation. |T^] One sees that 
the present theory is the non-Abelian generalization of that physics; in particular m, given above, 
is recognized as the Debye screening length, which remains gauge invariant in the non-Abelian 
context. 

It is especially interesting to emphasize that Silin did not use quantum field theory in his 
derivation; rather he employed classical transport theory. Nevertheless, his final result coincides 
with what here has been developed from a quantal framework. This raises the possibility that the 
non-Abelian Kubo equation can also be derived classically, and indeed such a derivation has been 
given, as mentioned above. 

We now pause in our discussion of solutions to the non-Abelian Kubo equation in order to 
describe its classical derivation. 

Transport theory is formulated in terms of a single-particle distribution function / on phase 
space. In the Abelian case, / depends on position {x^} and momentum {p^} of the particle. For 
the non-Abelian theory it is necessary to take into account the fact that the particle’s non-Abelian 
charge {Q“} also is a dynamical variable: satisfies an evolution equation (see below) and is an 

element of phase space. Therefore, the non-Abelian distribution function depends on {x^}, {p^} 
and {Q^}, and in the collisionless approximation obeys the transport equation -^f = 0, i.e. 

df df dp^ df dQ'^ ^ 

dx>^ dr ^p^^ dr dQ°^ dr 
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The derivatives of the phase-space variables are given by the Wong equations, for a particle with 
mass 


dx^ 

dr 

dr 

dQ^ 

dr 


pfi 


jptiu dxy 

fabcdx^ Abn 


c 


In order to close the system we need an equation for . In a microscopic description (with a single 
particle) one would have = / dxQ'^Y Y d‘^(x — x(t)) and consistency would require 

covariant conservation of the current; this is ensured provided satisfies the equation given above. 
In our macroscopic, statistical derivation, the current is given in terms of the distribution function, 
so the system of equations closes with 

= I dpdQQ^p‘^fix,p,Q) 

The collisionless transport equation, with the equations of motion inserted, is called the Boltzmann 
equation. The closed system formed by the latter supplemented with the Yang-Mills equation 
is known as the non-Abelian Vlasov equations. To make progress, this highly non-linear set of 
equations is approximated by expanding around the equilibrium form for /, 


free 

f 

■'boson 

fermion 


OC 




This comprises the Vlasov approximation, and readily leads to the non-Abelian Kubo equation. |]_ 

One may say that the non-Abelian theory is the minimal elaboration of the Abelian case needed 
to preserve non-Abelian gauge invariance. The fact that classical reasoning can reproduce quantal 
results is presumably related to the fact that the quantum theory makes use of the (resummed) 
1-loop approximation, which is frequently reorganized as an essentially classical effect. Evidently, 
the quantum fluctuations included in the hard thermal loops coincide with thermal fluctuations. 

Returning now to our summary of the solutions to the non-Abelian Kubo equation that have 
been obtained thus far, we mention first that the static problem may be solved completely. [ 111 
When the Ansatz is made that the vector potential is time independent, A± = A±(r), one may 
solve for a± to find a± = —A± and the induced current is explicitly computed as 


m 


■fj, 

.^induced 



0 ) 


This exhibits gauge-invariant electric screening with Debye mass m. One may also search for 
localized static solutions to the Kubo equation, but one finds only infinite energy solutions, carrying 
a point-magnetic monopole singularity at the origin. Thus there are no plasma solitons in high-T 
QCD.[1T| Specifically, one finds that upon selecting the radially symmetric solution, which decreases 
at large distances, there arises a magnetic monopole-like singularity at the origin. 

Much less is known concerning time-dependent solutions. Blaizot and Iancu|^] have made the 
Ansatz that the vector potentials depend only on the combination x ■ k, where k is an arbitrary 4- 
vector: A± = A±{x-k). Once again a± can be determined; one finds a± = and the induced 
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current is computable. For k 
time is present) one finds 



, where there is no space dependence (only a dependence on 


m 


•fl 

^ induced 





More complicated expressions hold with general k. The Kubo equation can be solved numerically; 
the resulting profile is a non-Abelian generalization of a plasma plane wave. 

The physics of all these solutions, as well as of other, still undiscovered ones, remains to be 
elucidated, and I invite any of you to join in this interesting task. 
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